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Abstract. We consider the behaviour of weak solutions of the unforced three- 
dimensional Navier-Stokes equations, under the assumption that the initial condition 
has finite energy (||u||^ = / but infinite enstrophy (||_Du||^ = Jjcurlup). We 
show that this has to be reflected in the solution for small times, so that in particular 
+00 as t ^ 0. We also give some limitations on this 'backwards blowup', 
and give an elementary proof that the upper box-counting dimension of the set of 
singular times can be no larger than one half. Although similar in flavour, this 
final result neither implies nor is implied by Scheffer's result that the 1/2-dimensional 
Hausdorff measure of the singular set is zero. 
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Introduction 

The existence or otherwise of singularities in the three-dimensional Navier-Stokes 
equations is acknowledged to be one of the outstanding open problems in the theory 
of partial differential equations, and so there has been much research concerned with 
the time before the formation of a possible singularity. In this paper we consider the 
behaviour of solutions after the occurrence of a singularity, or equivalently (but less 
hypothetically) the time evolution of a weak solution with a suitably irregular initial 
condition. 

In order to describe our results in more detail we need some simple terminology and 
notation. Throughout the paper we consider the evolution of an incompressible fluid in 
a periodic box Q = [0, L]^ in the absence of external forces: the evolution of the velocity 
field u{x) is modelled by the three-dimensional Navier-Stokes equations 

— - uAu + (u ■ W)u + Vp = V • w = 0. 
ot 

For convenience sake we also assume that the total momentum is zero, u{x) dx — 0, 
a situation that is preserved under the evolution of u. 
We denote the kinetic energy of u{x) by 

\\u\\'^ := / |M(a:)|^dx 
Jn 

and the 'enstrophy' of u{x), which is the square integral of the vorticity u!{x) — curl u{x), 

by 

:= / |curlw(x)|^dx. 
Jn 

We denote the space of all divergence-free flow fields with finite kinetic energy by H, 
and the space of all such fiow fields with finite enstrophy by V. (More strictly, if is 
the space of all C°° Q-periodic functions that are divergence free and have zero average 
over Q, then H is the closure of Y in the norm, and V is the closure of Y in the 
norm.) 

Leray showed in 1934 that given an initial condition with finite energy, there is a 
weak solution of the equations that exists for all time and whose kinetic energy remains 
finite. Indeed, we have 



u{t)f+ f \\Du{s)fds<\\uof 
Jo 



'0 

for all t > 0, and the solution is weakly continuous into H, i.e. u{t) u{s) as i — > s. 
However, the uniqueness of such a solution is not guaranteed. 

On the other hand, while it is know that if the enstrophy ||DM(i)|p remains finite 
then the solution is unique, an initial condition uq with ||Dmo|P < cxo is only known 
to give rise to a solution whose enstrophy remains finite for a limited time. [For an 
elementary discussion of these issues see Robinson (2006); for a more advanced but 
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accessible treatment see Doering & Gibbon; and for a more technical mathematical 
discussion see Constantin & Foias (1988) or Temam (1977).] 

The possibility of singular behaviour in the Navier-Stokes equations was first 
pointed out by Leray (1934) who suspected the existence of self-similar singular 
solutions. He also showed (using a simple argument which we will review later) that 
any solution that blows up at a time t = T must be bounded below by 

\\ou(t)r > 

Still little is known of the regularity (or otherwise) of weak solutions of the Navier- 
Stokes equations, but several papers have imposed some restrictions on the hypothetical 
blow up. The most famous is probably Scheffer's result stating that the Hausdorff 
measure of turbulent times in the Navier-Stokes equations is less than one half. More 
sophisticated is the result of Caffarelli, Kohn, & Nirenberg (1982), which treats the 
singular set in space-time (i.e. as a subset of i7 x [0, oo)). Scheffer's result can be 
obtained as a consequence of the result of Caffarelli et al., see Section 7 in Robinson 
(2006). 

Here we give a simple proof of a result reminiscent of Scheffer's: we show that the 
upper box-counting dimension of the singular set is no greater than 1/2. This result 
neither implies, nor is implied by, Scheffer's. 



1. Solutions with irregular initial data 

Given initial data Uq G V, it is well known that the solution remains regular (i.e. \\Du{ 
remains finite) at least for some short time interval. Indeed, standard estimates can be 
used to show that for t > s 

\\Duit)r< , (1) 

- ^\-c{t-s)\\Du{sW 
while the right-hand side remains finite, i.e. while 

c{t-s)\\Du{s)f<\. (2) 

In particular, if s = then this shows that the solution remains regular at for t G [0, T), 
where T = c~^||Duo||^^. 

Rearranging ([1]) to obtain a lower bound for s < t gives 



Pn(s)f >^^i^^=. (3) 

- ^\^c{t-s)\\Du{tW 

While one might be concerned as to the range of validity of this bound (given the 
conditions on t and s in ([2]), it does in fact hold for all s <t. Indeed, if for some s <t 
one has 



\Du{s\\?< 



Du 



^yi + c{t - s)\\Du 
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then the right-hand side of ([T]) remains finite on (s, t] and one can infer that < 
||Du(t)||^, a contradiction. 

From ([3]) one can easily deduce (as did Leray) that if a solution does blow up as 
t — > T, i.e. if [[/^^(t)!! ^ oo as t — T, then for s < T one must have the lower bound 

\\Du{s)f> ^ for all < s < T. (4) 

V c(T - s) 

Given mq G V, we know that even if there is a loss of regularity at some time T, 
there is a weak solution that exists for all positive times. It is therefore reasonable 
to ask about the behaviour of this solution after the putative blowup time. Even the 
equations are in fact 'regular' - i.e. if any initial condition in V gives rise to a regular 
solutions that is valid for all time - one can still ask the same question about solutions 
whose initial condition is contained in H and not in V. 

Our first simple result shows that such solutions retain a trace of this 'irregularity', 
i.e. it is not possible that ||-Dmo|| = oo but is uniformly bounded on an interval 

of the form (0, e) for some e > 0. 

Lemma 1.1. Suppose that uq E H\V , and that u{t) is a weak solution with m(0) = Uq. 
Then \\Du(t)\\ —>■ +oo as t ^ 0. 

Proof. Suppose not. Then for some M > we have a sequence t„ with < M. 

Since u is weakly continuous into H, we know that u{t„) uq in H. However, since 
||-DM(tri)|| < M it follows from the Alaoglu Theorem that u{tn) Vq in V, for some 
Vo with ll-Dfoll < M. Since weak convergence in V implies weak convergence in H, 
and weak limits are unique, we must in fact have Mq = Vq, and so H-DmoII ^ M, a 
contradiction. □ 



2. Decay of solutions 

Given the result of this lemma, it becomes natural to look for possible upper or lower 
bounds on the rate of decay of the V norm for solutions with irregular initial data, and 
we now turn to some results in this direction. 

All of our results will be consequences of the lower bound in (I3l), which is more 
convenient to work with when written in the form 

\\Du{s)f>—, ^ . (5) 

- ^\\Dum-' + c{t-s) 

We note here that this lower bound remains valid even if u is not regular on the interval 
[s,t], since ||Z^M(r)|| = +oo for r G (s, t) implies that 

\\Du{s)r > > ^ 
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2.1. Solutions m L4(0,r;1/) 

An immediate consequence of the inequality in ([5]) is an upper bound on the V norm of 
solutions known to be in L^(0, T; V^). In order for 



\\Du{s)rAs 



to be finite, one would expect blows up more slowly than t~^ as t ^ 0. In 

particular, it is natural to imagine that remains bounded as t — > 0, and this 

is shown rigorously in the following result: 

Proposition 2.1. Suppose that u G L^(0,T;y). Then there exists a constant c 
depending on ||M||L4(o,T;y) such that 

t'l^\\Du{t)\;' < c. 



(6) 



Proof. From ([3]) we have 



\\Du{s)r> 



\\Dum 



l + c{t-s)\\Du 



and so 



WuisWds > 



— \og(l + c(t - s)\\Du 
c 

log{l + ct\\Du{t)\\^) 



J 



from which (l6l) follows. 



□ 



It is natural to conjecture that this decay also holds for all weak solutions (this 
is the case for the two-dimensional Navier-Stokes equations), but in general all the 
information we currently have available is that u G L'^{0, T; V), and so a more natural 
conjecture along similar lines is that < c. However, a similar approach, but 

starting from ([5]) yields much weaker results. Indeed, if = +oo with t„ | 

then one can only obtain 

oo 

a/2 



■n+lj 



< CXO, 



n=l 



an observation originally due to Leray. 



2.2. Box-counting dimension and decay of solutions 

Given a set X, its (upper) box-counting dimension is defined as follows. Let N{X, e) 
denote the largest number of disjoint balls of radius e with centres in X. Then 



df{X) 



lim sup 

e-*0 



\ogN{X,e) 
-loge 
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i.e. essentially N{X,e) ~ g-'^/C-^). [Usually one would take N{X,e) to be the smallest 
number of balls of radius e whose union covers X (see Definitions 3.1 in Falconer, 1990), 
but the two definitions are equivalent and that given here is more convenient for what 
follows.] 

We now give an indication of how the box-counting dimension occurs naturally 
when considering rates of decay of solutions. Suppose that 

pn(t)f>ct-". (7) 



Then, since we know that \\Du{s)\\'^ ds < ||moIP) it must be the case that a < 1. 
However, if (I7j) holds then in particular it follows that for t„ = (n/c)~^^" we have 
1 1 -Dm IP > n. The fractal dimension of this sequence is equal to q;/(1 + a), and 
since we should have a < 1, it follows that the dimension of this sequence should be 
less than one half. 

Rather than proving a strong estimate such as that in ([7j), instead we prove that if 
we have a sequence > n, then the fractal dimension of can be no larger 

than one half. 

Proposition 2.2. Suppose that 

\\Du{tr,)f > n. 

for some sequence tn > 0. Then df{{tn}) < 1/2. 
In the proof we write ^ = {tn}. 

Proof. Suppose that df{^) = d > 1/2. Then for some 6 with 1/2 < 6 < d there exists 
a sequence ej such that 

AT, :=iV(^,e,)>e-^ 

where N{^,e) denotes the maximal number of disjoint balls centred on tj of radius e. 
Let the centres of these balls be 1 < i < Nj. Since these balls are disjoint, we have 



Du{s)fds >J2 \\Duis)fds > 



i=i '"i <=j 

Now, using the bound in ([5]) it follows that 



\\Du{s)fds. 



\\Du{s)fds > 



t-e 



1 *"i 



ds 



t-e 



1 

2^ 
1 

2^ 



^\\Du{t)\\-'^ + c{t-s 
-^\\Dum-^ + c{t-s) 

^/\\Du{t)\\-^ + c€-\\Du 



t 

t-e 



-2 



Noting that VX2Tce - X is a decreasing function of X we therefore have 



/' 

^0 



\Du{s)\?ds > 



Ni 

2c ^ 



n 
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/ \\Du{s)rds>^Y.yn-^ + ce, 

"^0 n=l 

In order to proceed we note that 

X<v^/2 ^ VX2 + ce - X > 
It follows, since iVj > ej^ and 5 > 1/2 we have, for j large enough 



2(cej)-l/2 



^ n=2{c£j)-i/2 

Since ^ as j — oo and 5 > 1/2 it follows that /^^ ||-Dm(s)||^ ds = +oo, a 
contradiction. □ 

Our discussion above suggests that in fact one should expect to be able to prove 
that df{3^) < 1/2 provided one assumes in addition that tn i (an assumption not 
necessary for the proof of the Proposition). This is further supported by the following 
observation. If one has 

p?z(n-i)f > n, (9) 

i.e. tn = in the statement of Proposition 12.21 then although df{^) = 1/2 one can 
show that is impossible. Indeed, for all s G (tn+i^^n) we have 



vi+(i„-s)iiBK(i„)r ,/rT 



n{n+l) 

So 

/ \\Du{s)r>y / , , =>y^ — = +oo, 
Jo -^n(n+i)^/r7zzr-^y2(n + i) 



n{n+l) 

which is a contradiction. 



3. The set of singular times 

A very similar, but slightly simpler argument, shows that if S is the set of 'singular 
times' of a weak solution, i.e. 

S = {t > : u{t) ^V} = {t>0: \\Du{t)\\ = oo}, 

then < 1/2. 
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Corollary 3.1. The upper box-counting dimension of the set S of singular times satisfies 
< 1/2. 

Proof. Suppose that df^E) = d > 1/2. Then as in the proof of Proposition 12. 2[ for some 
5 with 1/2 < 5 < d there exists a sequence ej such that Nj = N{T,,ej) > ej^. If 
there centres of these balls are tn, 1 < n < Nj, then since the balls are disjoint 

/ \\Du{s)fds>Y, \\Du{s)f ds > J2 \\Du{s)fds. 

Since ||DM(t„)|| = +oo, the lower bound in (jll) holds, and so 

1 

Since the right-hand side tends to infinity as j oo, we must have dfiY?) < 1/2. □ 

Scheffer's well-known result that the 1/2-dimensional Hausdorff measure of the 
set of singular times is zero, je^/^{T.) = 0, is based on the same ingredients as our 
result, but neither is a consequence of the other. Indeed, in general the box-counting 
dimension gives an upper bound for the Hausdorff dimension, but from this observation 
the consequence of our result is only that ^^/^(S) < +oo. However, note that while 
any countable set K has Hausdorff dimension zero (so in particular Jif^lK) = for all 
s > 0), the box-counting dimension has the property that df{X) = df{X), so it is easy 
to find a countable subset of the line with df{X) = 1. Thus our result df(T,) < 1/2 does 
indeed serve to limit further the set of singular times. 
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